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Abstract 

A proposal is made for a mathematically unambiguous treatment 
of evolution in the presence of closed timelike curves. In constrast 
to other proposals for handling the naively nonunitary evolution that 
is often present in such situations, this proposal is causal, linear in 
the initial density matrix and preserves probability. It provides a 
physically reasonable interpretation of invertible nonunitary evolution 
by redefining the final Hilbert space so that the evolution is unitary or 
equivalently by removing the nonunitary part of the evolution operator 
using a polar decomposition. 
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Evolution in spacetimes containing closed timelike curves has drawn much 
attention [1-8] recently following the work of Morris, Thorne and Yurtsever 
on the possibility of time machines [|]]. Such a spacetime is not globally hy- 
perbolic and therefore does not admit a foliation by spacelike hyper surf aces. 
Standard quantum field theory techniques cannot be applied. Nevertheless 
it has been proven that free field evolution on examples of such backgrounds 
is consistent and unitary while in the case of interacting fields, the evolution 
is apparently nonunitary[5-7]. The concept of nonunitary evolution is suf- 
ficiently disturbing that one might be tempted to dismiss these results out 
of hand. This would be a mistake as a closer study is repaid by a clearer 
understanding of unitary evolution. 

There have been previous considerations of nonunitary evolution and pro- 
posals for how to interpret it. Jacobson argued that causality implies uni- 
tarity, and that there are ambiguities in computing expectation values for 
observables localized in regions spacelike separated from a region of nonuni- 
tary evolution ||. The possibility was raised that a nonunitary Heisenberg 
formulation could be causal and unambiguous. 

An unambiguous but acausal path integral interpretation based on the 
decoherence functional formalism was proposed by Friedman, Papastama- 
tiou and Simon|Q. This was elaborated thoroughly by HartleQ. Both find 
that the presence of closed timelike curves in the future acausally affect the 
probabilities for observations in the present. In addition, the decoherence 
functional is nonlinear in the initial density matrix, and this represents a 
nonlinear (but unobservable) modification of standard quantum mechanics. 

In this paper, after reviewing the discussions of Jacobson and Hartle, a 
proposal is made for an unambiguous and causal treatment of (invertible) 
nonunitary evolution which is linear in the initial density matrix. This pro- 
posal is given both in the language of states and as a decoherence functional. 
The resulting description in fact shows that the evolution is unitary when 
appropriately handled. 

The essence of the proposal follows from the observation that in [5-7] the 
evolution X is identified as nonunitary because XX* ^ 1, but this is not 
sufficient to conclude nonunitarity in the sense of loss of probability |TIJ. A 



unitary operator U is defined[ll] to be a linear operator from one Hilbert 



space Ti. a onto a second Hb such that 

m\u4>) Hb = {^\<t>) Ha , (i) 

for all -0, G 7i a . Because all Hilbert spaces with a countable basis are 
isomorphic |TT| , it is conventional to write inner products formally without 
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stating their measure densities explicitly. Additionally, in standard quantum 
mechanical examples, the Hilbert space obtained after evolution is the same 
as the initial one and the inner products are the same. This combination of 
facts makes it easy to overlook that the final inner product and Hilbert space 
must be specified when defining the adjoint. (X-0|X0) Wa and {Xif:\Xip) Hb 
are not the same unless Ti a = TCb- Thus, X may not be unitary (X^X ^ 1) 
when the final Hilbert space is taken to be the same as the initial one, but 
nevertheless X is unitary in the sense of (|I]) for a suitably chosen final Hilbert 
space. The proposal for handling (invertible) naively nonunitary X is to 
work with the final Hilbert space where X is unitary. This is equivalent 
to remaining in the original Hilbert space 7i a but removing the nonunitary 
part of X through a polar decomposition X = RUx, where Ux is unitary 
UxUx = 1 and R 2 = XX^ a is self-adjoint in 7i a . Then, Ux is the manifestly 
unitary evolution operator in 7i a . This will be discussed in detail below. 

One's initial concern with nonunitary evolution is that probability must 
not be conserved, and that therefore one cannot compute expectation values. 
A proposal to get around this is to simply renormalize expectation values by 
the (changing) non-unit norm of the state in which the expectation is taken 

This gives a well-defined expectation value even when the norm of is 
not unity. The difficulty with this proposal, as shown by Jacobson||, is that 
it can give ambiguous results. 

Consider a measurement of P a made in a region A well-separated in a 
spacelike direction from a region of nonunitary evolution, for instance a region 
containing closed timelike curves (CTCs). The expectation value of P a can 
be computed on either a spacelike hypersurface cr_ through A that passes 
before the CTC region or on a spacelike hypersurface a + through A that 
passes after the CTC region. The wavefunction \ifj) is different on these two 
hypersurfaces and is related by the nonunitary evolution X N through the 
CTC region 

|V>) + = x„N>)_ (3) 

where the subscript labels the hypersurface on which the wavefunction is 
evaluated. The expectation value of P a on <r_ is then 

(P) - -MW>- (4) 
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while on cr+, it is 

Since the observable P a is measured in the same region on both slices, it has 
not changed under evolution, that is, Xj^ P a X^ = P a - This implies that 

= Mx n pm 

{ °' + -Wx'nXnM- 1 

which is not equal to (P a )- unless X^ N X N = 1. Thus, for this proposal, 
nonunitary evolution leads to ambiguous predictions for expectation values of 
observables even in regions spacelike separated from the region of nonunitary 
evolution. 

A different proposal M using a decoherence functional in Gell-Mann and 



Hartle's generalized quantum mechanics fll2|, [13[ succeeds in avoiding this 
ambiguity. The decoherence functional measures the interference between 
histories characterized by strings of projections. The projections serve to 
divide the set of fine-grained histories into coarser-grained classes. The deco- 
herence functional can be computed as a closed-time-path functional integral 
over histories which start from an initial density matrix, pass forward through 
one sequence of projections up to a final surface, and then back through a 
different set. With no final condition, one traces over a final surface. In 
the Heisenberg picture of the unitary case, projection operators on a space- 
like slice (Tk can be expressed in terms of the unitary evolution U(ak,o~ a ) of 
projection operators defined on an initial slice <7j 

P^k) = U- 1 (* k ,* i )PtoU(<r k ,* i ). (7) 

An overall unitary operator evolving from the initial surface cr, to the final 
surface tjf cancels by cyclicity of the trace and the decoherence functional is 
left as a trace over the initial hypersurface 

D(a', a) = Tr(p^K) • ■ • P$MPiP™M ■ ■ ■ P^M) • (8) 

In the non-unitary case discussed above, there is a region bounded by 
a + and o~_ containing a non- unitary evolution X^. This evolution must be 
included in the decoherence functional: Xn is inserted in the sequence of 
projections where it takes place. Projections before <r_ are expressed in the 
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Heisenberg picture in terms of the initial hypersurface <7j, and those after a + 
in terms of projections on the final hypersurface <7y 



U- l (a k ,a f )P^U(a k ,a f ), a k > a+. 



(9) 



To make the transition between the initial and final surfaces, one uses the 
full evolution 



Additionally, the decoherence functional must be renormalized to account for 
the nonunitary evolution of the initial density matrix in the absence of coarse- 
graining. The proposed decoherence functional in the case of one projection 
made before the CTC region is|| 



For the case of one projection made after the CTC region, it would be 



Now, if the projection is made in a region spacelike separated from the CTC 
region, then as above P a commutes with X and one sees that D_(a',a) = 
D + (a',a). There is no ambiguity It has disappeared because both compu- 
tations are aware of the full future history, and in particular of the presence 
of the CTC region. 

This proposal has a curious property however. The result of computing 
the decoherence functional is different in the presence of X than it is in its 
absence. If we are living in the past of a CTC region, we might naively com- 
pute with a decoherence function which did not include X. Our predictions 
would be different than they would be if we included a nonunitary X. If 
the correct physical description involved X, the difference between physical 
results and our no-X predictions would in principle be susceptible to exper- 
imental detection. With this decoherence functional, the future existence of 
a CTC region acausally affects the accuracy of predictions today. This is a 
strange result, but very important if true. 



X = U(a f ,a + )X N U(a_,a t ). 



(10) 




(11) 



D 




(12) 
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A second unusual feature of this decoherence functional (and indeed the 
earlier proposal) is that it is nonlinear in the initial density matrix. That 
is, because the trace of the initial density matrix appears in the denomi- 
nator, one cannot add different initial density matrices. This is therefore a 
nonlinear modification of conventional quantum mechanics. Arguably in a 
closed-system formulation of the wavefunction of the universe, we would not 
be sensitive to this in any event as we only see linearity in the superpositions 
of subsystems, each of which are part of the full initial density matrix. 

There is however an alternative. A new proposal leads to both unambigu- 
ous results and is causal and linear. This proposal is based on the observation 
that a linear invertible transformation X, such that XX^ a ^ 1, is not a uni- 
tary transformation from the Hilbert space 7i a to itself, but it is a unitary 



transformation to a new Hilbert space 7i&||10||. The states obtained by uni- 
tary evolution with X are elements of Tib, and they must be transformed 
back to 7i a if one is to compare with the physical states defined there. After 
this is done, one finds that these evolved states are simply those obtained 
by evolving with the unitary part Ux of a polar decomposition of X = RUx 
[where R = (XX^) 1 ^ 2 is self-adjoint in 7i a ]- Unitarity is restored by remov- 
ing the nonunitary part of X. This proposal can be implemented both in the 
language of states and expectation values and in the decoherence functional. 

The manipulations needed to elaborate the proposal can all be performed 
symbolically, but for clarity and familiarity, the case of quantum mechanics is 
treated explicitly first. In the language of states, one has to be careful about 
the form of the inner product when it is computed on different hypersurfaces 
after evolution. The inner product on an initial hypersurface <7j for states in 
the Hilbert space TC a is 



O*a)(0>#> 1 ') = / dx<f>*fi a r(j (13) 

= (l)(0,^a|^,«), 

where the (possibly operator-valued) measure density in the inner product 
is pre-subscripted parenthetically. (The measure density is put in front to 
emphasize that it is more naturally associated with the bra-vectors than the 
kets.) Let \ip, f) = X\ip,i) be the state on the final hypersurface a/ obtained 
by evolving an initial state on Oi with X — the tilde indicates that this is a 
state in the Hilbert space Tib- Requiring that the value of the inner product 
be preserved when evaluated on the hypersurface 07, one has 

(Ma)(^MlVM) = / cfe(A~ 1 0)*/i a X~ 1 '0 
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f 

= I dx^fiaX'^-X- 1 ^ (14) 
= / dx4>*^= M (<p,f\^,f), 

The induced measure density of the unitarily equivalent inner product is 

H„ = (X-y^aX- 1 = n a (XX^)-\ (15) 

The notation fi indicates that the adjoint is computed in the trivial measure 
density fi — 1. Note that this adjoint is related to the adjoint in the measure 
density \i a by 

X^fi a = ti a X u . (16) 

One sees that to preserve the inner product the measure density changes dur- 
ing the evolution, acquiring an inhomogeneous factor (II^)- 1 measuring 
the nonunitarity of X. 

A side remark may prevent potential confusion. Consider the (nonrela- 
tivistic) quantum mechanics of a particle propagating in a time-dependent 
curved background. The measure density /i a is the square-root of the de- 
terminant of the metric and therefore is time-dependent. Nevertheless it is 
specified once for the whole spacetime and the form of the inner product 
never changes; on each hypersurface where one computes the inner product, 
one evaluates fi a at that instant of time. Here, the transformation to /U& 
reflecting the nonunitarity of X is on top of any time-dependence \x a may 
already have. It is a genuine change to a different Hilbert space. 

Having found the Hilbert space in which X is unitary, we must face 
the fact that it is not the Hilbert space we naively expected. Our usual 
association between mathematical and physical quantities is made with the 
presumption of a particular Hilbert space. We know what physics means in 
the usual Hilbert space, and we would be surprised if it made a difference to 
be told that there had been "nonunitary" evolution in our past and that the 
Hilbert space at present is not the one we naively thought! Fortunately, we 
needn't worry: the choice of measure density/Hilbert space is on the same 
par as a choice of coordinates; you can freely change it when talking about 
physical quantities. We just need to find the transformation from /x b to [i a . 

Let R = {XX^y/ 2 . Then R = R u is self-adjoint in H a . Returning to 
the computation in Eq. (0), one sees that defining 

\ijJ} = R- 1 \^f}, (17) 
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resets the measure to /j, a , 



/ dxcf)* /j, a R 2 ip — dxtp* Haip . (18) 
f f 

But this means that the evolved state in 7i a is 

\^J) = R- 1 X\i;, i ) = U x \ij,i). (19) 

It is easy to check that U x is unitary in H a (U X U% = Uj?U x = 1)! This, X = 
RU X , is the polar decomposition of X in TC a . The restoration of unitarity is 
explicit. 

Under the transformation QT7| ) from Tib to 7i a on the final surface, the 
self-adjoint observable P a = P^ in Hb becomes 

P a = R~ l P a R, (20) 

which is self-adjoint in Ti a . With this in mind, one can turn to consider "non- 
unitary" evolution in the Heisenberg picture. The naive Heisenberg evolution 
of the operator P a is X _1 P a X. Unfortunately this cannot be correct because 
X^P a X is not self-adjoint in 7i a (that is, on the initial hypersurface). This 
should not be surprising. 

Heisenberg evolution works by pulling back operators from the final sur- 
face to the initial surface. If P a is self-adjoint in Tib on the final surface, then 
it is easy to confirm that its pull-back is self-adjoint on the initial surface in 

x- 1 P a x = x- 1 p£>x = X^^P^^X (21) 

= {x- l P a x)^. 

Using Eq. (|20|), this implies the Heisenberg evolution of P a is 

X- l P a X = U^P a U x . (22) 

Heisenberg evolution is accomplished using the unitary part of the polar 
decomposition of X. 

Sandwiching this Heisenberg operator between initial states, one sees that 
one has what one wants-the expectation of P a in 7i a on the final surface, 

M(il>,i\UZ l P a Ux\il>,i) = (^MflPaWJ). (23) 
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From another perspective this result is natural: in the Heisenberg picture, 
one does not have the freedom to change the Hilbert space, so to restore 
unitarity one must change one's identification of the evolved observable. 

One might be concerned that the nonlinear operation involved in remov- 
ing the nonunitary part would prevent one from composing X with a further 
unitary operator, say Uy. Let Z = UyX be such a composition, then 

u z = (zz u y 1/2 z 

= (UyR 2 Uy X )- 1/2 U y RU x (24) 
= U Y U X . 

Here one can use UyR^Uy 1 = (Uy RUy 1 ) 2 before taking the square root, but 
it is generally true that fiUyRUy 1 ) = Uyf(R)Uy l . Unitary composition 
works as one expects and there is no change to ordinary quantum theory in 
regions with unitary evolution. 

If the expectation value of an operator P a is computed on the initial 
hypersurface <7j, one finds 

Ui*)( P a) = (l)(^,i\^aPa\^,i)- (25) 

On the final hypersurface 07, the "push-forward" of P a is XPaX^ 1 which is 
assumed here to remain P a . Its expectation value on 0$ is then 

{ , b) (XP a X- 1 ) = M (P a ) (26) 
= w&flXPaX-^f) 

= {1) (4>,i\n a P a \4>,i). 

This agrees with the result computed on cr.j and there is no ambiguity. As 
well, the state does not appear in the denominator and expectation values 
take their usual linear quantum mechanics form. Finally, the result does not 
know of the presence or absence of X and is not acausally affected by CTC 
regions. 

Turn now to quantum field theory. To quantize a field, one begins in a 
non-interacting initial region and solves the (initial) free wave operator Lj 
for its modes, say L^Uj = 0. The modes are normalized in an inner product 
of the form fllcf ). The wave operator is self-adjoint in the measure density 
of the inner product. Using a scalar field as an example, the quantum field 
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is expanded in terms of the modes with creation and annihilation operator 
coefficients 

'I' »», • ( 27 ) 

j 

The vacuum is defined as the state annihilated by all the cij, and the n- 
particle states are formed by applying a product of creation operators to the 
vacuum. 

The next step is to compute new mode function for the (final) free wave 
operator in a non-interacting final region. In the original measure density fi a , 
this operator is Lf, the usual free wave operator evaluated in the final region. 
It has modes v^. If a nonunitary evolution X takes place and the measure 
density in the final region has changed (0), then Lf is not the self-adjoint 
free wave operator in /aj, but instead Lf = RLfR' 1 is. The mode functions 
Vk of Lf are related to those of Lf by 

v k = (28) 

This corresponds to Eq. (|Tj]) above. 

Take the non-interacting case first. To compute the S-matrix, the initial 
modes Uj are evolved into the final region by X where they may then be 
expanded in terms of the complete set of modes Vk (and t>£). Bogoliubov 
coefficients are calculated in the usual way using the inner product with 
measure density ^b- The S-matrix is necessarily unitary because the inner 
products of states are preserved by the evolution and because the basis of final 
mode functions is complete. Alternatively, if one chooses to work always in 
the original familiar Hilbert space, one simply uses the final mode functions 
Vk, evolves the initial modes with Ux = R~ 1 X and computes Bogoliubov 
coefficients using the inner product with \x a . 

In the interacting case, attention shifts away from evolution of the mode 
functions to the detailed correspondence between initial and final n-particle 
states. The evolution operator X can be formally expressed as 

X= E a nm\n,f) {fla) (m,i\. (29) 

{n},{m} 

The n is short-hand for the detailed structure of the n-particle states associ- 
ated with the Vk modes, and m represents m-particle states associated with 
the Uj modes. If X is not unitary in 7i a , then 

XX U =i?Vl, (30) 
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where R = Ftf a is self-adjoint in 7i a . Multiplying this equation on the right 
and left by one finds it defines a unitary transition matrix Ux = R~ 1 X. 
This matrix is the S-matrix for the interacting theory and is unitary by 
construction. It is non-trivial if Ux is not the identity. Unitarity has again 
been restored by removing the nonunitary part of the evolution operator. 
For completeness, I remark that it follows from Eqs. (O) and ( p0|) that 

X h X = XX U = R 2 ^ 1. (31) 

Thus, R is self-adjoint in Tib as well as H a . Multiplying by R^ 1 on the left 
and the right, one finds that Ux = XR^ 1 is unitary in Tib- This explains the 
role of the second polar decomposition X = UxR- it provides the unitary 
operator for evolution which remains in Tib- 

This story can be retold in the decoherence functional formulation. There, 
unitarity is restored by appropriate definition of the trace. The trace in the 
decoherence functional is taken over a final hypersurface. One must make 
sure that it is the trace which is physically equivalent to one defined on an 
initial surface, say before the CTC region. This can be done by requiring that 
the trace of pi is preserved under the evolution. On an initial hypersurface 
o"j, the measure density is p a and the trace is 

Tl (Ma)(pi) = Tr ( i)(p^ a ). (32) 

On the final hypersurface cr/, the measure density changes to pb and one 
takes the trace of the evolved density matrix XpiX' b 

Tr {tlb) (X Pi X^) = Tr {1) (X Pi p b X^) (33) 
= Tr (1) (X PiX^fjL b ). 

The traces of pi on oi and <jf are equal if 

p b = (X-^paX- 1 = pa(XX^)- 1 = p a (X^X)-\ (34) 

as above. It is important to emphasize that care must be taken with the 
use of the different types of adjoint as one changes the hypersurface/Hilbert 
space where one is taking the trace and as one adjusts whether the measure 
density is implicit or explicit. 

I remark that the trace of p\ is also preserved 

Tr (Ma)(P* 2 ) = Tr {1) {pip a pip a ) (35) 
= Ti {l) {X Pl X^pbXp t X^p b ) 
= Tr M ((Xp, t X^) 2 ). 
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Thus, this proposal does not provide an interpretation for evolutions which 
are not invertible in which a pure state evolves to a mixed state. Such an 
evolution would not be described by XpiX^, but the point is emphasized so 
that there should be no mistake. 

The decoherence functional for one projection made on the hypersurfaces 
<j\ and o"2 before and after the CTC region is given by 



D(a',P';a,P) = (36) 
Tr (w )(p5 ) (^ 2 )XP^Vi)A^ 1)t6 (^)^^ 2)t6 (^; 



[ (i) 

The projection operators are in the Heisenberg picture as defined in Eq. 
(|9p. This means that P^ (tfi) is self-adjoint in TC a on the initial hypersurface 
while P%f > {o2) is self-adjoint in Ti.b on the final hypersurface. In the first form, 
the trace is evaluated on the final surface in Tib. The second form makes the 



dependence on the measure density explicit. This satisfies Hartle's axioms |L3 
for a decoherence functional. 

By using Eq. fl3~4|), the decoherence functional can be expressed in terms 
of a trace in Ti, a on the initial hypersurface, 

D(a',(3';a } f3) = (37) 
= Tr {1) (p£\a 2 )XP! 3 }\a l )p t P^^ 

The effect of the change in measure has been to pull back the projection 
operator P^ 2 )(cr 2 ) to the initial hypersurface by the adjoint action of X^ 1 . 
Since (X^ 1 )^ 1 p a X = fi a R~ 2 , the decoherence functional can be reexpressed 
as a trace in 7i a on the final hypersurface 

D(a',P';a,P) = (38) 

= Tr (1) (i^V2)XP£Vi^^^^ 

There is a factor R~ 2 present here that is not in the naive decoherence 
functional, cf. the numerator of Eq. ([12|) . A more subtle but equally im- 
portant difference is that the projection operator after the CTC region is 
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Pa 2 (°2) which is self-adjoint in Jib, not 7iJ Using Eq. ( |20"| ) to transform to 
the projection operator P^ 2 )(a 2 ) which is self-adjoint in H a , the expression 
for the decoherence functional is found 

D(a',(3';a : (3) = Tr (Ma) {PS\a 2 )U x P^\a l )p l P^\a l )U^ P^{a 2 )) . 

(39) 

This clearly shows that the mechanism for unitarity restoration is to replace 
the nonunitary evolution X by its unitary part Ux- 

Finally, the trace in the decoherence functional can be shifted to the initial 
hypersurface by inserting UxUx 1 = 1 at the end of the trace. Using cyclicity, 
one sees that the decoherence functional is built from the Heisenberg evolved 
projection operators (|22| ) 

D(a',P';a,P) = (40) 

= Tr (Ma) ([/^ 

Working with Eq. (p7|), one can consider a few cases to discover the 
dependence of the decoherence functional on X. Suppose first that there is 
only a projection before the CTC region so that P^(a 2 ) and P a i(<72) are 
not present. Then, the X factors cancel and there is no dependence on X. 
A CTC region in the future does not acausally affect prediction made in the 
present. Alternatively, suppose that P^(<J 2 ) and Pffi (cr 2 ) are in a region 
spacelike separated from the CTC region so that they commute with X. 
Again the X factors cancel and the CTC region has no acausal effect. Only 
if P^\g2) does not commute with X will there be an effect. This presumably 
can happen if P^(a 2 ) is in the causal future of the CTC region. 

One can also compute probabilities for alternatives in the CTC region 
by inserting projection operators between periods of nonunitary evolution. 
These alternatives must of course be members of an exhaustive and exclusive 
set. Because there is not a spacelike hypersurface which intersects the CTC 
region, it is not clear how to describe such a computation in terms of states 
at the instant projection is made. The computation can be made in the 
Heisenberg picture by pulling the projection operator back to the initial 
hypersurface and computing the expectation value there, as in Eq. fl3~T|). 

It should be mentioned that one cannot necessarily detect CTC regions 
in one's past because experimentally we only have access to an effective 
initial density matrix. Our effective initial density matrix may have been 
renormalized by X from the true initial density matrix, but we would not 
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be able to detect this without knowing the initial density matrix. On the 
other hand, if one has a theory for the initial condition of the universe, as in 
the Hartle-Hawking no-boundary proposal, predictions based on this initial 
condition could differ depending on the presence or absence of CTC regions 
and this might be detectable. Otherwise, the best means of detecting a CTC 
region is to measure a change in amplitudes as one evolves past it-in this 
way, one could compare results before and after and see a possible difference 
due to the noncommutativity with X. Sadly, there seems little likelihood of 
an experimental test of this kind. 

For completeness, one can consider a time-symmetric decoherence func- 
tional in which a final state density matrix is also imposed. The objective 
here is to verify the time symmetry of the decoherence functional in the 
presence of nonunitary evolution. The subtlety in this is again taking care 
with the hypersurface on which the trace is evaluated. To be completely 
explicit, consider a situation in which an evolution X from an initial surface 
a a to an intermediate surface o\ is followed by a projection F^ 1 ^ and then by 
subsequent evolution Y to a final surface o\>. (Note that in contrast to the 
discussion above, the projection operator is in the Schrodinger picture and 
the intermediate evolutions are included explicitly as X and Y.) Neither 
X nor Y will be assumed to be unitary. To be clear about the definition 
of adjoints, let fi a be the (coordinate- valued) measure density on the initial 
surface a a . The adjoint f a on a a is defined by \i a X^ a = X^p, a and the adjoint 
ffe is defined similarly on the final surface using the measure density fib- The 
projection is self-adjoint = Pa ^ m on the intermediate surface where it 



is defined, using the induced measure density there |14 



On the final hypersurface the decoherence functional is 

Tr M (p f YP^X Pi X^P^Y^ 

D(a';a) = ^ , (41) 

Tr M (p f YX Pi XW>>) 

Tr (1) (pf^YP^XpiX^P^Y^^ 



Tr^pf^YXpiX^Y^^ 

The final density matrix acquires a factor of the measure density when one 
introduces the measure explicitly. This is easily confirmed with a calculation 
in terms of states where a trace over two density matrices leads to a product of 
two inner products and hence requires two measure density factors. Using the 
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transformation formula for the measure density, = (YX) lj!l fx a (YX) x , 
the denominator is seen to be symmetric 

Tr {1) [p f ^YX Pi X^Y^^ = Tr^^X-'Y-'pfiY-y^X- 1 )^^. 

(42) 

The time-symmetry of the numerator is slightly more difficult to establish 
due to subtlety with the adjoints. Two side results are needed. First, the 
inverse and adjoint (fi) operations commute 

(X-y 1 = (X h )-\ (43) 

which follows by taking the adjoint of XX~ X = 1. Next, self-adjointness 
of P^ on the intermediate surface implies self-adjointness of the projection 
pulled back to the initial surface X~ 1 P^X when the measure density trans- 
forms appropriately for nonunitary X, cf. Eq. (j2ip . Note however that 
pCO ^ pWt« if x is not unitary 

In the numerator of (|4"ID, consider the factor ^YP^pX. This can be 
reexpressed with \x a on the right on transforming the measure and using the 



computation fl2ll) 



^YXix^p^x) = ((yx)- 1 )^^ 1 ^?^ 1 ^- 1 )^ 

= (y-^pW^X- 1 )^. (44) 

Similarly, one has 

XtipWt^ti^ = ^x-ipWy-i. (45) 

These can be combined to show that the numerator of the final equality in 
©is 

Trw^/^-^^F^^X-^tv^X^pWy- 1 ) = (46) 



This is the time-symmetric form of the numerator, and one concludes that 
the decoherence functional with initial and final density matrices is time- 
symmetric. 

To summarize, it has been shown that evolution by an invertible "nonuni- 
tary" operator X with XX^ a = R 2 ^ 1 can be made unitary by transforming 
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the measure density on the final hypersurface appropriately, Eq. ([TJ^). Equiv- 
alently, one can stay in the original Hilbert space by removing the nonunitary 
part of X to obtain the unitary evolution operator Ux = R~ l X. Both proce- 
dures lead to unambiguous results for the expectation values of observables. 
The results are causal and do not involve nonlinear modifications of quantum 
mechanics. This proposal for handling nonunitary evolution can be stated in 
the form of a decoherence functional in the language of Gell-Mann and Har- 
tle generalized quantum mechanics. The decoherence functional with initial 
and final density matrices imposed as conditions is time-symmetric. 

The conclusion is that one can consistently make quantum mechanical 
computations in spacetimes containing closed timelike curves without loss of 
probability or causality if evolution is described by an invertible nonunitary 
operator. A region of closed timelike curves could only be detected if one 
interacted with something in its causal future. 
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